function [c,ceq] = mycon(x)
% function [c,ceq] = mycon(x)
% constraint function for penalty_bc_shock.m in fmincon, restricts the
% vector to have length 1

c = [];         % Compute nonlinear inequalities at x.
ceq = x'*x-1;   % Compute nonlinear equalities at x.
